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ABSTRACT 

Generalized Hankel transforms are useful in analyzing the effect of circularly 
symmetric optical systems on arbitrary inputs. Some examples of such systems are 
complex laser resonators and space telescopes. Three methods for performing Hankel 
transforms with optical or digital processors are described. The first method is 
applicable when the input data is available in Cartesian (x-y) format and uses the 
close connection between generalized Hankel transform and the two-dimensional 
Fourier transform in Cartesian coordinates. The second method is useful when the 
input data is in polar (r - 0) format and uses change of variables to perform the 
nth order Hankel transform as a correlation integral. The third method utilizes 
the von Neumann addition theorem for Bessel functions to extract the Hankel 
coefficients from a correlation between the radial part of the input and a Bessel 
function. Initial experimental results obtained for optical implementation of the 
first two methods are presented. 


INTRODUCTION 

The analysis of complex optical systems is greatly facilitated by two- 
dimensional Fourier transform techniques. The effect of an optical system on 
arbitrary inputs is easily described by a transfer function in the Fourier domain. 
Generalized Hankel transforms are similarly useful when dealing with a circularly 
symmetric (or axisymmetric) system for arbitrary inputs.^ This situation is 
encountered in performing mode analysis on the output of a slightly misaligned 
laser resonator as well as in aligning space telescopes. An optical method for 
performing mode analysis via generalized Hankel transform will have the unique 
advantage of preserving the phase of the wavefront to be analyzed. 

It is well known that when a two-dimensional function has circular symmetry 
(i.e. it depends only on the radial variable, r) , its Fourier transform is also 
circularly S 3 mmietric (i.e. it depends only on the radial variable, p) . It can be 
shown that in such a case the Fourier transform is equivalent to the order 
Hankel transform . ^ 
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where { } indicates two-dimensional Fourier transformation. 
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where ? (p) is the 0 order Hankel transform and J^(pr) is the 0^^ order Bessel 
function of the first kind. Thus in dealing with circularly S 3 nnmetric systems, the 
Qth order Hankel transform (which is a one-dimensional operation) can be used 
instead of the two-dimensional Fourier transform if the inputs also are circularly 
symmetric. For arbitrary inputs, this technique can be extended by using the genera-- 
lized Hankel transform, which expands on the nth order transform (F„(p) with kernel 
(p) in Eq. (D). ^ " 

The generalized Hankel transform, t>e defined for an arbitrary 

function, f(r,0), as follows: 
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The generalized Hankel transform thus involves a Fourier series expansion in 0 , 
followed by an n^^ order Hankel transform (with n^^ order Bessel function, J^(pr), 
as the transformation kernel) on the n^^ coefficient of expansion, f^(r). This 
generalized Hankel transform is useful in analyzing systems with circular symmetry 
but when the input is not circularly symmetric. 

In the following sections we will describe three techniques for easy 
implementation of the generalized Hankel transform with optical or digital 
processors. These three techniques are each applicable in different circum- 
stances. Initial experimental results on the optical implementation of the 
first two methods will be presented. In conclusion, we will detail the course 
of future work in this area. 

OPTICAL IMPLEMENTATION OF GENERALIZED HANKEL TRANSFORM 

In optical processors a two-dimensional Fourier transform with respect to the 
Cartesian coordinates (x,y) is performed very easily with the help of a simple 
spherical lens.^ The equivalence of two-dimensional Fourier transform and Qth 
order Hankel transform for circularly S 3 mimetric functions was established in the 
Introduction. A similar connection exists between two-dimensional Fourier trans- 
form with respect to Cartesian coordinates (x,y) and n^^ order Hankel transform 
of the radial part if the function is of a* form f(r) ed^^ . 
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where 

2 2 ^^ -1 
p = (u +v ) , (() = tan (v/u) . 


and by definition 
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This connection arises out of the integral representation of (pr) . ^ Using this 
result in the definition of generalized Hankel transform given in Eq* (2) the 
following relation is obtained between the two-dimensional Fourier transform and 
the generalized Hankel transform: 


oo 



t tl 

Thus the generalized Hankel transform F^^^Cp) is equivalent to the n^ 
coefficient of Fourier series expansion of ^^(p,(|>) in variable (|) . Using 
this result, an optical system shown in Fig, 1 is designed to perform 
generalized Hankel transform on an arbitrary input that is available in 
Cartesian format (f(x,y)). The spherical lens performs a two-dimensional 
Fourier transform on f(x,y) generating ^^(u,v) • A suitably designed computer 
generated hologram then performs the Cartesian-to-polar (u,v->p,(})) coordinate 
transformation on^^(u,v) generating ^^(p ,({)), ^ This is followed by a cylindri- 
cal lens which takes a one-dimensional Fourier transform with respect to variable 
(|) generating the desired output, 

Initial optical experiments established the connection between two- 
dimensional Fourier transform and n^^ order Hankel transform for a function 
of the form f(r)eJ^^, The schematic diagram of the optical system is shown 
in Fig. 2. A computer generated hologram was used to encode eJ^® dependence 
of the input. The output was detected by a TV camera, which measures the light 
intensity in the Fourier plane of the input. The input used in these experi- 
ments had an r dependence given by 6(r-a), thus corresponding to a thin ring 
of radius "a" in the Cartesian (x,y) plane. The n^^ order Hankel transform of 
6(r-a) is aJ^(ap), making the output easily understandable. Figure 3 shows the 
results corresponding to 0^^ order Hankel transform (i.e, eJ^^ =1), In Fig, 3(b) 
the function | J^Cap) P is plotted, which is then compared to a line scan through 
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the origin (Fig. 3(c)) of the output shown in Fig. 3(a). 


Very good qualitative 
Figures 


agreement between the theoretical and experimental results is obtained 
4(a-c) present the results for I.q order Hankel transform for the same f(r), 
but here the 0 dependence is e^ . Again good qualitative agreement is 
seen between theory and experiment. A line scan through the origin of the output 
corresponding to 2^^ order Hankel transform of 6(r-a) (with 0 dependence e3^®) is 
shown in Fig. 5. The zero at the origin was broader and the side lobes were seen 
to fall slower indicating that we indeed have as expected. 

OPTICAL IMPLEMENTATION OF N™ ORDER HANKEL TRANSFORM 


The previous method is applicable when the input is available in Cartesian 
(x,y) format, since it involves performing two*~dimensional Fourier transforms 
with respect to the Cartesian variables. If the input is polar (r,0) formatted, 
a more direct approach outlined in Eq. (3) has to be followed in obtaining 
generalized Hankel transform. The first part of the operation, which involves a 
Fourier series expansion in variable 0 , is easily performed optically using a 
cylindrical lens. The calculation of n^^ order Hankel transform of the n^“ 
coefficient of expansion is less straightforward since it corresponds to a 
space-variant operation. So the main aim of the next two methods is to perform 
the space-variant operation of n^^ order Hankel transform optically. 

a. Method Using Change of Variables 

One standard procedure used in converting a space-variant operation into a 
shift- invariant operation is to employ appropriate change of variables. In the 
case of n^b order Hankel transforms the following procedure was described by 
Siegman for implementing the space-variant operation as a correlation integral 
on a digital processor.^ From the definition 
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Fn(y) = ^ f(x) J^(xH-y) dx 
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where F^(y) = pF^(p), f(x) = rf(r), and J^(x+y) = arp J^(rp). The algorithm, 
therefore, consists of first linearly weighting the input f(r) and perform- 
ing (logarithmic related) coordinate transform. This distorted input 
is then correlated with a similarly weighted and coordinate-transformed n^h 
order Bessel function to give the desired Hankel transform as well as linearly 
weighted and coordinate distorted forms. In any physical system the correla- 
tion integral will be performed over a finite interval, giving rise to 
truncation errors. Also if the input is sampled in the x-domain, the sampling 
rate should be adequate to represent the function accurately in x-domain. 
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These factors and others are discussed at length, especially for digital 
implementation, in Refs. 4 and 5. 

Since the operations of coordinate transformation and correlation can be 
performed by an optical processor, the optical system outlined in Rig. 6 can 
calculate n^^ order Hankel transforms. Computer generated holograms are used 
to perform r->x coordinate transformation as well as to encode the Fourier plane 
filter with impulse response j^Cx). The second dimension of the optical processor 
can be used to perform different order Hankel transforms on different inputs, thus 
achieving multichannel operation. 

In the initial experiment, both the input and the Fourier plane filter were 
encoded by computer generated holograms. The linear weighting and the coordinate 
transformation of the input was performed by the digital computer before doing the 
holographic encoding. The computer generated holograms used the Lee~Burckhardt 
technique^ and contained 128 pixels. The optical system is depicted in Fig. 7. 

This system was then used to perform 0^^ order Hankel transform on two different 
inputs, fi(r) = PiJo(Pl^) f2(^) P2'^o(P2^)* result of 0^^ order Hankel 

transform on fqCr) and 12^^) should be 6 (p-p 2 ) and 6 (p-p 2 ) respectively. The 
results of the computer simulation of this algorithm are depicted in Fig. 8 . The 
finite width of the peak and the sidelobes are due to the finite limits of intergra- 
tion. To facilitate easy comparison with the experimental results, the order 

Hankel transform with linear weighting and coordinate distortion (i.e. jF^Cy))^) 

was plotted versus y instead of F (p) versus p. Figure 9 shows the experimental 
results obtained. The optical output was detected by a 1024-element Reticon 
linear photodiode array. The shift in the peak position and the difference in the 
peak heights (due to the linear weight) are evident, indicating good qualitative 
agreement with the computer simulation results. 

b. Method Using Neumann Addition Theorem 

This method investigates an approach based on the special properties of Bessel 
functions.^ If the input f(r) is correlated with Jj^^(pr) then the m^^ order Hankel 
coefficient F^(p) is obtained at the origin of the correlation plane for the 

particular value of p encoded in the Bessel function kernel. 
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The Neumann addition theorem for Bessel functions states that 


( 6 ) 


J (p(r+r’)) = J (pr’) J (pr) (7) 

m m^n n 

n=-oo 

Substituting for J (p(r+r’)) from Eq. (7) into Eq. ( 6 ) we get 
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Thus it is seen that the correlation plane contains an infinite sum of Bessel 
functions of different order weighted by Hankel transform coefficients of 
different order. So in principle it is possible to extract Hankel transform 
coefficients of different order out of a single one-dimensional correlation 
operation. 

The other dimension of optical system can be used to perform correlations 
with J^(pr) with different values of p to obtain complete Hankel transform. This 
approach is currently being evaluated further to determine the situations in which 
this method will be suitable. 


FUTURE WORK 

In this paper we have presented brief outlines of three different approaches 
for performing generalized Hankel transforms by optical or digital processors. 

We also presented initial experimental results on two of the three approaches. 

The future work will concentrate on carrying out the coordinate transformation 
required by the first two methods optically via computer generated holograms. A 
more quantitative analysis of the performance of the optical processor will be 
carried out. These three methods will be compared with each other with respect 
to space-bandwidth requirements as well as pre- and post-processing requirements. 
Finally, generalized Hankel transforms will be applied to specific problems, such 
as analysis of complex laser resonators, and optical systems will be developed for 
those problems. 
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Figure 1. - Schematic diagram of optical processor for performing generalized 
Hankel transform via a two-dimensional Fourier transform with respect to 
Cartesian coordinates (x,y) . 




■F^ {f(r) ei""} = F„(p) e”'"'* 

Figure 2* - Schematic diagram of experimental setup for establishing the 

th 

relation between n order Hankel transform of f(r) and two-dimensional 
Fourier transform of f(r) e"^^^ with respect to Cartesian coordinates (x>y)* 
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(a) Photograph of the output* 



(b) Plot of |j (ap) I versus p, which is 
‘o' 

the theoretically expected result. 


(c) Line scan through the origin of the 
pattern in figure 3 (a) , giving 

If (p) 1 versus p. 
o 


Figure 3 . 


Experimental 


results for 0 

F (p) of 6(r-a) 
o 


order Hankel transform 
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Figure 5. - Line scan through the origin of the two-dimensional 


Fourier transform of 6(r - a) e 
6 (r - a) . 
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Figure &• ~ Schematic diagram of an optical processor for performing 
order Hankel transform on f(r) employing change of variables. 














(b) Input J^(P 2 r). 


Figure 9. - Oscilloscope traces of output of the optical processor 

j . 1 — 

performing the 0 order Hankel transform. Traces correspond 

, /N I 2 

to |F (y) I versus y for the two inputs. 
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